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SUMMARY 


Empirical- analytic methods are presented for estimating pressure and 
thermal loads induced by trailing-edge eleven deflections on adjacent surfaces 
of high-speed aircraft. The methods are applicable for turbulent-boundary-layer 
flews. The extent of the shock-wave — boundary-layer interaction region may be 
estimated adequately for all planar surfaces, but not for curved surfaces. How- 
ever, knowing the extent of the three-dimensional inter action- flow region, pres- 
sure and heating distributions are adequately predictable. Heating magnification 
factors can be estimated using the methods presented herein once the pressure 
distributions are known. 


INTRODUCTION 

Interactions between the shock wave and boundary layer can appreciably 
affect the pressure and heating distributions on high-speed aircraft surfaces 
(refs. 1 to 16). The changing loads can affect aircraft stability and alter 
control effectiveness . Also, shock-wave — boundary- layer interaction flows can 
cause very high aerodynamic heating to occur at local surface regions. Any and 
all of these inter act ion- flow facets can compromise an aircraft design. 

Except for simple two-dimensional (2-D) shapes, theoretical methods are not 
available for predicting either the changes in the surface pressure and thermal 
loads, or even the extent of the interaction. Not only is the three-dimensional 
(3-D) inter act ion- flow problem of more pr act i cal design importance, but it is 
further exacerbated by both the increased magnification of 3-D interaction- flow 
effects, and the lack of theoretical methods which would enable designers to 
predict these severe loads. 

A photograph of a hypersonic research airplane model is shown in f igur e 1. 
Three-dimensional interaction- flow regions of paramount importance are expected 
in the engine inlet area and where shock waves from the fuselage , wings, fins, 
and elevens impinge on other surfaces of the aircraft. These inter action- flow 
regions can result in very large local-pressure and thermal peaks. Unfortu- 
nately, cctnplex 3-D interaction- flow regions are not amenable to purely theoret- 
ical analyses (refs. 1 to 20) . 

Recourse must be made to empirical methods for the design of high-speed 
aircraft. A model of a hypersonic research airplane is shown in figure 1 and 
described in reference 6 . Many organizations have done a great deal of work 
toward obtaining an understanding of the driving mechanisms of 3-D shock-wave— 
boundary-layer interaction flows (refs. 1 to 20) . There have been many experi- 
mental investigations of these complex flows for both total-aircraft conf igura- 
tions (ref. 6) and also for fundamental model configurations representative of 
component parts of hypersonic aircraft. Ex am pi es of the latter are steps or 
fins mounted on flat plates or bodies of revolution (refs. 1 to 5 and 1 4) . 



The particular 3-D interaction flow addressed herein is that caused by 
deflecting trailing-edge elevens. The shock wave generated by the eleven can 
impinge on an aft fuselage surface and the shock-wave — boundary- layer interac- 
tion can result in a much larger region of disturbed surface pressures and heat- 
transfer distributions than anticipated using inviscid-flow analyses. Sketches 
of a typical hypersonic research airplane and a wing-elevon model, representa- 
tive of the aft portion of the airplane, are shown in figure 2. Detailed data 
obtained for shock-wave interactions on the wing-elevon model (ref. 2) were 
very helpful in developing methods for estimating surface pressure and thermal 
loads in interaction- flow regions. 

A related problem that has received much attention is that caused by fin 
generated shock waves incident to a boundary layer on an adjacent planar surface 
(refs . 8 and 12). The prime difference being the existence of a boundary layer 
upstream of the elevon hinge line. The similarity in the interaction flows 
prompts the use of both types of data in developing methods for estimating 
3-D interaction- flow effects. 

Many empirical methods that had been postulated for predicting surface 
pressure and thermal loads in 3-D interaction-flow regions were examined and 
compared with experimental results from other investigations. Portions of the 
more promising methods were modified slightly to improve agreement with an 
extended base of experimental data. The engineering methods presented herein 
were developed using extensive data from many sources. They are the best 
methods available for calculating pressure and thermal loads caused by shock 
waves, induced by elevon deflections, interacting with turbulent boundary 
layers on adjacent surfaces of hypersonic aircraft. 


SYMBOLS 

a location of virtual origin along ^-axis (eg. (7) and fig. 11 ) , cm 

b semiminor axis defining hyperbola shown in figure 11, cm 

b/a slope of asymptote to separation location (eg. (6) and fig. 11) 

c empirical proportionality factor for peak heating (eg. (1 6) ) 

HL hinge line 

2 

h heat-transf er coefficient, W/m -K 

L distance f ran plate leading edge to elevon hinge line (eg. (17)), cm 

£ separation length upstream of hinge line (fig. 7) 

M Mach number 

n exponent in expression for peak pressure magnif ication (eg. (1 3) ) 

p pressure. Pa 


2 



R unit Reynolds number 

X streamwise distance measured downstream frcsn where elevon-generated 

shock wave intersects wing surface (fig. 11), cm 

x' streamwise distance measured downstream from wing leading edge along 

surface of wing and elevon, cm 

z distance measured on end-plate surface upward from wing surface, cm 

Az distance from shock wave to onset of pressure rise (eg. (8)) 

6 boundary-layer thickness, cm 

e elevon deflection angle, deg 

6 elevon shock-wave angle, deg 

\ angle frcm vertex of shock and x-axis to ray delineating maximum 

pressure or heating (eg. (10)), deg 

shock-oriented, rectangular coordinate system (eg. (5) and fig. 11), 

Subscripts: 

incip incipient 

max maximum value 

peak peak value 

plat plateau value 

sep separation 

u undisturbed 

0 undisturbed flow value at hinge-line location 

1 local flow conditions over wing and end-plate surfaces undisturbed 

by elevens 

2 two-dimensional, obligue-shock value 

<» free-stream flow conditions 

TYPES OF INTERACTION FLOWS 

Two types of 3-D interaction flows induced by trailing-edge elevon deflec 
tions are sketched in figure 3. For smaller elevon deflection angles, the 
elevon-generated shock emanates from the hinge line (fig. 3(a)). For larger 


cm 
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deflection angles, the flow separates from the wing surface upstream of the 
hinge line, and a shock emanates from the separation line (fig. 3(b)). In both 
cases, the extent of the disturbance on the adjacent end-plate surface is far 
larger than would be anticipated using inviscid-flow analyses. The extent of 
the disturbed flow is shown in figure 3 by oil-accumulation lines. 


WING AND ELEVON SURFACES 
Pressure Distributions 

Deflection- angle effects .- The local-flow compression angle is of first- 
order importance in determining the extent of 2-D turbulent boundary- layer 
separation ahead of elevens (flaps or ramps) and the resulting pressure dis- 
tribution (refs. 21 to 29) . The pressure rise required to cause 2-D turbulent 
boundary-layer separation may be expressed as a weak function of the local 
running-length Reynolds number. For a wide range of supersonic/hypersonic 
test conditions, the onset of separation ahead of elevens of sufficient chord 
length that the flow reattaches on the eleven surface (Chord length » 6 q) 
occurs for compression angles between 1 5° and 20°, and is predictable using 
Korkegi's criteria (ref. 13) 


Pincip 

Pi 


1 + 0.3M-|^ 


0.1 7Mi 

L 


for Ml < 4.5 


for M-) >4.5 


( 1 ) 


More recent methods show excellent correlations between incipient separation 
and friction coefficient. The small range of angles, the weak dependence on 
Reynolds number, the unsteadiness of turbulent separation, and the existence 
of 3-D eddies in 2-D separated flows are factors that seem to make superfluous 
more sophisticated means for predicting the onset of 2-D turbulent boundary- 
layer separation than that in equation (1) (refs. 21 to 29). 

Example schlieren photographs of 2-D turbulent boundary-layer flows ahead 
of elevons are shown in figure 4 (ref. 2) . The free-stream Mach number is 6 
and the free-stream Reynolds number, based on the distance frem the unswept wing 
leading edge to the eleven hinge line (64.14 cm), is 18.5 million (ref. 2). The 
local Mach number over the wing surface upstream of the eleven is 5.9. The 
shock wave emanating from the eleven hinge line changes slightly in shape from 
convex at low eleven angles to concave at larger eleven angles. A small region 
of separated flow in the vicinity of the hinge line is discernible only for a 
30° eleven deflection. This is in agreement with the Korkegi criteria (ref. 1 3) 
shown in equation (1 ) . 

Sample pressure distributions, obtained in the Langley 20-Inch Mach 6 Tunnel 
(ref. 17), are compared with inviscid flow calculations in figure 5. Turbulent 
boundary layers do not separate significantly upstream of elevons for pressure 
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rises less than those calculable using equation (1) (i,e., no separation for 
G < 26° based on inviscid-flow analysis and the conditions shown in fig. 5). 

As would be expected, the pressure distributions on the elevon surfaces are sub- 
stantially different from those that would be calculated using inviscid-flow 
theory. For supersonic flows, the pressure on the flap remains less than the 
inviscid value for many boundary-layer thicknesses (refs. 21 to 29). For 
hypersonic flows however, there are generally "overshoots" in the flap-surface 
pressure distribution for sufficiently large deflection angles (i.e., angles 
large enough for separation) . The loads on the flaps may be estimated assuming 
that the pressure attains the inviscid value at a surface distance of five 
boundary-layer thicknesses downstream of the hinge line (see "Wing and Control 
Loads" section). 

Wing-sweep effects .- As long as the boundary layer is fully turbulent at 
all spanwise stations upstream of the hinge line, wing sweep has a negligible 
effect on the surface pressure distribution on the wing and elevon (ref. 1 7) . 
There is no significant spanwise change in the pressure distribution. For 
trailing-edge controls that cause separation ahead of the hinge line, the 
extent of separation is uniform in length upstream of the hinge line for fully 
turbulent boundary layers . It is nonuniform if the boundary layer on the wing 
surface is transitional or laminar on the outboard portion of a swept wing where 
the chordwise Reynolds number is smaller (refs. 30 and 31 ) . 

Center -body and wing-tip-fin effects .- Vertical streamwise planar surfaces, 
such as end plates and tip fins, have little influence on the pressure distribu- 
tion over the wing and elevon surfaces for attached flows (refs. 32 and 33) . 
However , if the elevon angle is large enough to cause separation, end plates 
prevent venting of the vortical separated flow and the extent of turbulent 
boundary-layer separation ahead of trailing-edge elevens is more nearly uniform 
(ref. 32) . For fully turbulent boundary-layer separation, the addition of end 
plates does not significantly extend the region of separated flow ahead of a 
trailing-edge elevon (compare figs . 83 and 84 in ref . 34) . However , if the 
boundary layer is laminar or transitional near the separation location, the addi- 
tion of end plates greatly enlarges the extent of separation, changes the effec- 
tive surface shape, and results in far different pressure distr ibutions than 
obtained without end plates (refs. 32 to 34) . 

A nonplanar center body (fuselage) has a negligible effect on the pressure 
distribution on the wing and elevon surfaces as long as there is no flow separa- 
tion (ref. 32) . However , even without separation, a fuselage causes an outflow 
of oil-flow streaks on the wing surface (ref. 32 and fig. 6(a)). When the elevon 
angle is large enough to cause turbulent boundary-layer separation, the fuselage 
shock and pressure fields cause large changes in the pressure distributions and 
extent of separation on the wing and elevon surfaces. As evidenced by the oil- 
flow photographs in figure 6, the maximum extent of separation ahead of the 
elevon is twice as large when the cylindrical body is attached to the wing as 
when the end plate is attached to the wing. As would be expected, oil accumu- 
lates at the juncture of the silicone rubber insert and the stainless steel wing, 
even though this juncture is thoroughly sealed. Nevertheless, the juncture is 
well upstream of the interaction region. 
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A planar end plate on the inboard portion of the wing, or a wing-tip fin, 
influences the separated flow slightly in its immediate vicinity (ref. 32 and 
fig. 6). However, a nonplanar fuselage causes a remarkable change in the 
separated- flow region, as evidenced by the oil-flow photographs shown in fig- 
ure 6. The spanwise location of the maximum extent of separation occurs where 
the fuselage shock (calculated using the method described in ref. 35) intersects 
the eleven hinge line (fig. 6(c)). Corresponding to the larger extent of 
separation on the wing surface, reattachment occurs farther aft on the eleven 
surface. Although the spanwise location of this maximum extent of separation 
is predictable (occurring where the fuselage shock intersects the hinge line) , 
there is no method for predicting the chordwise extent of separation. Recourse 
must be made to experimental results. 


Wing and Control Loads 

Wing loads .- Loads on the wing surface are unaffected unless the trailing- 
edge control causes the boundary layer to separate from the wing surface. 
Pressures are increased in the region of turbulent boundary-layer separation 
ahead of a control. A great amount of data have been reviewed in order to 
obtain methods for estimating the increased load on the wing surface when turbu- 
lent boundary- layer separation occurs upstream of the hinge line. Empirically, 
a dividing streamline angle of approximately 1 5° is observed for turbulent 
separation at Mach numbers from 2 to 1 0 and local running-length Reynolds num- 
bers frem 5 million to 50 million. Based on a review of extensive data (e.g., 
refs. 15, 16, and 30), the following equation was developed for estimating 
increased pressure levels on wings in regions of turbulent boundary-layer 
separation upstream of elevons: 


— = 1 .8 + O.llM-j ^ (2) 

Pi 


This empirical equation fits available data frexn many sources, and within the 
scatter of the data (20 percent) the pressure rise is nearly independent of 
local running-length Reynolds number for fully turbulent boundary-layer 
separation. 

The extent of 2-D turbulent boundary-layer separation depends primarily 
on the flap deflection angle and the chord of the trailing-edge flap. The 
separated flow usually reattaches either on the flap surface or at the trailing 
edge of the flap. In the case where the flow reattaches on the flap surface, 
the extent of separation scales with the undisturbed boundary-layer thickness 
6 and the unit Reynolds number. In the case of reattachment at the flap trail- 
ing edge, the upstream extent of separation scales with the flap chord. 

Ranges of upstream lengths of separation, nondimensionalized with respect 
to undisturbed turbulent boundary-layer thicknesses, are shown in figure 7 for 
a wide range of flow conditions (Mach numbers from 2 to 1 0 and local running- 
length Reynolds numbers from 5 million to 50 million) for cases where reattach- 
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merit occurs on the eleven surface (refs. 23 to 29). Generally, the normalized 
extent of separation increases with boundary-layer thickness and eleven deflec- 
tion angle, but decreases with increasing Mach number and Reynolds number. 

Indeed, more detailed correlations of separation lengths with Reynolds number 
or local skin-friction coefficients have been presented in references 23 to 29. 
These more complex correlations should be consulted for more specific estimates 
of separation lengths for cases where reattachment occurs on the flap surface. 

Whether the separated flow reattaches on the eleven surface or at the trail- 
ing edge of the eleven may be approximated with confidence by using inviscid- 
flow theory for sonic flow downstream of the eleven hinge line. When the calcu- 
lated inviscid flow over the eleven surface is subsonic, reattachment occurs at 
the eleven trailing edge. In this case, the extent of separation may be esti- 
mated using a dividing streamline angle of 1 5° with reattachment at the eleven 
trailing edge. For smaller eleven angles, where the single-shock inviscid flow 
over the eleven surface would be supersonic, the extent of separation over the 
wing surface may be estimated roughly using figure 7. The extent of separation, 
along with the increased pressure on the wing surface in the separated flow 
region (eg. (2)), may be used to estimate the additional force on the wing 
surface. 

Eleven loads .- When the eleven angle is sufficient to cause separation, the 
pressure on the eleven surface usually remains less than that which would be 
calculated inviscidly. Even in cases where there is no flow separation upstream 
of the hinge line, loads on the eleven surface vary from those that would be 
calculated using inviscid-flow methods. In general, the pressure on the eleven 
surface does not attain the inviscid pressure rise until 5 to 1 0 boundary-layer 
thicknesses downstream of the hinge line (refs. 23 to 29). However, for some 
instances of very high-speed fluid flow there is an "overshoot" in the sur- 
face pressure distribution (fig. 5 and ref. 36). In these cases the reflection 
waves caused by the shock wave emanating from the separation location should 
be considered. Otherwise, it is sufficient to assume that the pressure rise 
is linear to the inviscid value at a distance of 5 to 1 0 boundary-layer thick- 
nesses downstream of the hinge line. 

As evident in figure 5, the control force of an elevon is less than would 
be anticipated using inviscid-flow theory. This loss of control effectiveness 
is aggravated when the flow separates from the wing surface, increasing the load 
on the wing surface upstream of the hinge line (fig. 5) . 


Heat-Transfer Distributions 

Deflection-angle effects .- Heat-transfer coefficients on the wing and 
undeflected elevon surfaces may be calculated accurately (to within approxi- 
mately 10 percent) using any one of several methods (refs. 37 to 40). The heat- 
ing on most of the wing surface is not changed as long as the elevon angle is 
insufficient to cause the boundary layer to separate (ref. 32). However, in 
detailed experiments with solid cast models of a hypersonic research airplane, 
Lawing (ref. 6) observed very high heating along elevon hinge lines even in the 
absence of significant flow separation. 
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Heat-transfer -coefficient distributions on the eleven surfaces may be esti- 
mated assuming a virtual origin for the turbulent boundary layer at the hinge 
line and using the measured or calculated pressure distribution. This simple 
method, however, does not account for the thickness of the upstream boundary 
layer and can lead to large discrepancies between calculated and measured heat- 
transfer coefficients on the elevon surface (ref. 32) . 

Heat- transfer coefficients increase in regions of turbulent boundary-layer 
separation (refs, 14, 15, 16, 28, and 32). The increased heating corresponds 
closely to the increased pressures on the wing surface for turbulent boundary- 
layer separation; i.e., the extent of increased heating is the same as the 
extent of increased pressures. Heat- transfer coefficients are largest near 
the forward position of the turbulent separated- flow region (fig. 8) . The mag- 
nitude of these heating rates may be best estimated by comparing the design 
data with data gleaned from similar designs and flow conditions. 


Maximum heat-transfer coefficients occur near reattachment of the separated 
flow on the elevon surface (fig. 8). For turbulent separated flows, the maximum 
heat-transfer coefficients at reattachment may be estimated using (refs. 18, 28, 
and 41 to 44) 


'max 



0.85 


(3) 



This equation is helpful in estimating the order of magnitude of the maximum 
heat-transfer coefficient in most cases. However, measured peak heat-transfer 
coefficients can differ substantially from those calculated (ref. 40) , as evi- 
denced by the specific case shown in figure 8. As attested by many investi- 
gators (e.g., refs. 3, 5, 14, 15, and 16), there is considerable scatter and 
much uncertainty in heat-transfer data, particularly in separated- flow regions. 

The heat-transf er coefficients indicated by the solid line in figure 8 
were calculated using the measured pressure distribution over the elevon surface 
(refs. 32 and 40) . As is the case for the distributions shown in figure 8, the 
actual location of maximum heating is frequently upstream of that calculated 
from the measured pressure distribution. Farther downstream of reattachment, 
however, the calculated heat-transfer coefficients agree fairly well with those 
measured exper imentally . 

Wing-sweep effects .- Lines of equal heat-transfer coefficients parallel the 
wing leading edge (refs. 32 and 45) . As long as there is no flow separation, 
the heat-transf er distribution on the wing surface remains unchanged whether or 
not the elevon is deflected, when a turbulent boundary layer separates, the mag- 
nitude of the heat-transfer coefficient is increased in the separated- flow 
region. The area of increased heating is the same as that of increased pressure 
in the separ ated- flow region. 
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For swept wings, the heat- transfer cxiefficient is larger on the outboard 
portion of a deflected eleven. It is believed that this effect is due to the 
thinner boundary layer on the outboard portion of the swept wing (refs. 32 
and 45) . Since there may be a large difference between predicted and experi- 
enced values, within the broad range of Mach and Reynolds numbers stated previ- 
ously, the most reliable expression currently available for predicting the 
maximum heat-transf er coefficient seems to be that given in equation (3) . 

Center-body and wing-tip-fin effects .- Fuselages and wing-tip fins have 
relatively little influence on the heat-transf er distribution on the wing sur- 
face as long as there is no shock impingement on the wing. Again, the region of 
increased heating is the same as the region of increased pressure. The largest 
heat-transf er coefficients occur just downstream of the onset of turbulent 
k» un da ry- layer separation fron the wing surface (refs. 32 and 45) . 


VERTICAL SURFACES 

A deflected eleven, or trailing-edge flap induces pressure and thermal loads 
on an adjacent surface (such as an end plate, aft fuselage section, or vertical 
fin) over a much larger area than predicted by inviscid-oblique-shock analysis 
(see fig. 9) . The 3-D interaction flow induced by a deflected eleven is scxnewhat 
similar to that induced by a fin mounted normal to a surface, the difference 
being the existence of a boundary layer on the wing surface upstream of the ele- 
ven, leading to thicker boundary layers on the eleven surface and the possibility 
of flow separation upstream of the eleven. Conprehensive discussions of the 
nature of 3-D flow interactions are presented in references 5, 7, 8, 10, 11, 12, 
and 2 0 . 


Interaction-Flow Region 

Incipient separation .- Hayes (ref. 8), Korkegi (ref. 13), and others have 
pointed out that 3-D flew separation occurs much more readily than 2-D flow 
separation. The turbulent boundary layer on a surface subjected to a skewed 
shock wave will separate for streamwise pressure rises of 


Pincip 

« T ,65 

Pi 


(4) 


This value (1 .65) is simply an average of the values presented by Hayes in 
reference 8 (1 .8) and by Kor kegi in reference 1 3 (1.5). The onset of 3-D flow 
separation is difficult to as certain and has little effect on the loads on the 
end- pi ate surface until the pressure rise exceeds a factor of two. 

Pressure and thermal- load characteristics .- References 8, 41 , and 42 have 
proposed analytical methods for estimating pressure and heating distributions 
induced on a planar surface by a vertical fin with a sharp leading edge . Except 
for the boundary layer on the wing surface upstream of the eleven, their geom- 
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etry is similar to that discussed herein, with the fin simulating the elevon and 
the flat plate simulating the vertical end plate. Their methods are based on 
experimental results and correlations of empirical data frcan many sources. 

Sketches of interaction- flow regions on the end- plate surface and the cor- 
responding pressure and heat- transfer-coefficient distributions, as depicted by 
Hayes (ref. 8), are shown in figure 10. Surface streamlines, as evidenced by 
oil-flow investigations, start to curve and the pressure and heat- transfer coef- 
ficient start to rise at the onset of the interaction region. At separation, 
the pressure approaches a plateau value whereas the heat-transfer coefficient 
attains a local maximum value. At the shock location, the pressure dips to a 
local minimum value, then increases to a peak value near the elevon surface as 
shown in figure 10 (refs. 8, 12, and 41). Hayes (ref. 8) and Scuderi (ref. 41) 
observed that peak heat- transfer coefficients occur quite close to where the 
pressure attains its peak value, several boundary-layer thicknesses downstream 
of the elevon hinge line (fin leading edge). 

Extent of separation and interaction-flow regions .- Hayes (ref. 8) presents 
a method for estimating the location of separation that is found to be reliable. 
A shock-oriented coordinate system (^,C) is used with a virtual origin upstream 
of the shock wave emanating from the separation location ahead of the elevon 
(fig. 11). The location of separation on the adjacent end plate is well approx- 
imated by a branch of the hyperbola indicated in figure 11. Using the shock- 
oriented coordinate system and data from many sources, Hayes found that the 
data followed a hyperbola 





(5) 


as a good approximation for the location of separation. 

The slope of the asymptote to the separation location, gleaned frcan the 
results of many experiments, was presented in figure 1 0 of reference 8. Using 
the suggested parameter (M] sin 6), a simple equation was fit to the data pre- 
sented by Hayes 


b 0.968 0.127 

- = 0.442 (Ml sin 6) - 1 .262 + + (6) 

a (M] sin 0 ) (Ml sin 6 ) ^ 


where 6 is the shock-wave angle and M-) is the Mach number of the local 
undisturbed flow. The location of the virtual origin a along the C-axis 
may be estimated using the following equation, which is simply another curve 
fit to the extensive data presented by Hayes (ref. 8); 
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a 


5 


(7) 


= 3 + 


(M-] sin 6 - 1 . 3) ^ ^ 


where the distance a is in centimeters (curve fit to fig. 11 of ref. 8). 

The location of the onset of the inter act ion- flow region depends on both 
the shock strength and the local-boundary-layer thickness. Scuderi (ref. 41) 
relates the outboard extent of the beginning of the interaction pressure rise 
to the location of the shock wave, the pressure rise (governed by the eleven 
deflection angle e), and the boundary-layer thickness 6 


Az X 

— = (0.0115G + 0.1)- + 0.14e (8) 

6 6 


where G is in degrees and x is the streamwise distance downstream of where 
the eleven generated shock wave intersects the wing surface. This expression 
is generally valid and simple to use. The distance along the z-axis frem the 
shock wave to the onset of the pressure rise is Az measured downstream of 
where the shock wave intersects the wing surface (figs. 10 and 11). Hayes 
(ref. 8) presents data relating the pressure rise from the onset of the dis- 
turbed region to the location of separation, with the extent of the distance 
between the lines labeled "Onset" and "Separation" in figure 10. The over- 
pressure in this region, which averages approximately 0»5(pp]^g^. - p^ ) , adds 
little to the induced side force on the end plate. It may be accurately approxi- 
mated using a linear pressure distribution froti p-| at the onset location to 
Pplat separation location. 

Separation locations calculated using Hayes' method (ref. 8) agree well 
with experimental observations (ref. 32) shown in figure 12. The rectangular 
coordinate systan used in equations (5) to (7) and in figure 11 is shock wave 
oriented. The transformation equations are 


X = 5 cos 6 - C sin 6 - a cos 0 
z = K sin 6 + C cos 0 - a sin 6 


(9) 


where 6 is the shock-wave angle and the x, z origin is where the shock wave 
intersects the wing surface (not necessarily at the hinge line). 


Locations of peak pressure and heating .- Peak pressure ratios and maximum 
heat-transf er-ooef f icient ratios occur near the eleven surface. Hayes (ref. 8) 

X 

notes that for values of - < 10, the ray of peak pressure ratio is closer 

6 
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to the elevon surface than the ray of maximum heat-transf er ratio. Scuderi 
(ref. 41) has correlated much data «hd presents 


X - £ = 0.28(6 - e) for peak pressure 
X - e = 0.24(6 ~ e) for peak heating 


( 10 ) 


where X is the angle from the vertex of the shock wave and x-axis to the 
ray delineating the peak pressure or heating. Evident from equation (10), 
both of these rays are close to the elevon surface. For =5.9 
and e = 20°, 0 = 28.6°. From equation (10), (X - e)peak pressure ~ 2.4° 
and (X - e)pea)^ heating =2.1 . These are indeed close to the elevon surface. 


Pressure Distributions 

Plateau and separation pressure rises .- Hayes (ref. 8) and Scuderi (ref. 41) 
observed that pressure rises in the interaction-flow region correlate well with 
(M-| sin 0). A simple polynomial curve was fit to the data presented by Hayes in 
figure 7 of reference 8. The plateau pressure may be estimated using 


Pplat 

= 0.41 + 0.91 (Ml sin 6) - 0.06 (Mi sin 6)2 (11) 

Pi 


The pressure rise at the separation location may be estimated using data pre- 
sented by Hayes (ref. 8) 


Psep = 0-73Ppiat 


( 12 ) 


Scuderi (ref. 41 ) observed an initial peak pressure rise somewhat higher in 
value than the plateau pressure. This local pressure peak adds little to the 
side force acting on the end plate. In the experiments for M] = 5.9, no ini- 
tial peak pressure higher than the plateau pressure was apparent (ref. 1 7) . 
Neither Hayes (ref. 8) nor Korkegi (ref. 12) reported an initial peak pressure 
prior to the plateau-pressure region. 

Shock and peak pressures in inner-flow region .- The inner-flow region 
extends from the elevon surface to the shock wave (refs. 8 and 1 2) . There is 
a dip in the pressure distribution at the shock wave and then the pressure rises 
to a peak value near the elevon surface (refs. 8, 12, and 17) (fig. 13). The 
minimum pressure at the shock is difficult to ascertain, and the sparse data are 
difficult to correlate. The dip in the pressure apparently depends on shock 
strength, distance downstream of the shock vertex, and boundary-layer thickness. 
Although there is no convenient expression for the minimum pressure at the shock 
wave, the extent of the region is quiet small and the pressure ratio is bounded 


12 



by the plateau-pressure ratio and unity. The effect of this small dip in pres- 
sure on the side force exerted on the end plate is negligible. 

The magnitude of the peak pressure ratio E>peak/Pl also depends on shock 
strength, distance downstream, and boundary-layer thickness. However, the value 
of this ratio is well documented. The peak pressure is estimable using 
reference 8 


ppQ 3 k 

= (M] sin G) (13) 

PI 


by simply curve fitting equations to Hayes' data 


r 1.1 (x/6) 

1.4 + 

for 

X < 386 


1 3.8 + (x/6) 



(1 4) 

2.4 

V 

for 

X > 386 



Scuderi (ref. 41) also observed that the peak pressure increases with both shock 
strength and increased values of x/6. His results agree with those calculated 
using equations (13) and (14) and also with the results for M-j = 5.9 (ref. 17). 

Finally, at the elevon surface the pressure approaches the oblique-shock 
value obtained from reference 46 (refs. 8, 12, and 41) 


P2 7 (M-| sin 6)^-1 
PI 6 


(1 5) 


Close to the elevon hinge line (x < 56), the peak pressure may be less than the 
pressure at the elevon surface. Further downstream, the peak pressure exceeds 
the pressure at the elevon surface. 

Sample comparison .- A pressure-ratio distribution on the end-plate surface 
along a vertical line (calculated using the analytic method described in eqs. (5) 
to (15)), is compared with experimental data in figure 13. The particular case 
plotted is for M-| = 5.9, e = 20°, x = 7.6 cm, and 6 = 0.86 cm. The value of 
the pressure ratio at the elevon surface is that given by using equation (15). 

The peak pressure ratio is obtained using equations (1 3) and (1 4) and the loca- 
tion of the peak pressure is obtained using equation (1 0) . An S-shaped 

curve is faired between these two points in figure 13. The pressure on the 
end- plate surface, going away from the elevon, then dips to a minimum value 
(0.5(pp]^^^- + Pi ) ) at the inviscid-shock location (which is calculated using 
ref. 46). The pressure then rises to the plateau value (eq. (11)) and is essen- 



tially constant outboard to the separation location (calculated using eqs. (5) 
to (7)). Another S-curve represents the pressure-ratio distribution from the 
location of separation to the undisturbed value at the onset of the disturbance, 
given by equation (8) . Considering the complexity of the 3-D interaction-flow 
problem and the simplicity of the proposed analytical method, the agreement is 
quite adequate. 

Wing-sweep and tip-fin effects .- As long as the boundary layer is turbulent 
upstream of the elevon, wing sweep has no effect on the pressure on the end-plate 
surface. Tip fins at the outboard edge of the elevon usually have no effect on 
the end-plate-surface pressure distribution (refs. 1 7 and 32) . For very-low- 
aspect-ratio elevons, relatively large tip fins may result in more extensive 
separation ahead of elevons. This would shift forward the location of the 
interaction-flow region on the end-plate surface. 


Force on End Plate 

Elevons influence a much larger area of the end-plate surface than would be 
predicted using simple 2-D inviscid-flow analysis. However, the average pressure 
rise in this region is considerably less than that downstream of a 2-D oblique 
shock wave. In the experiments for = 6, the total side force induced on the 
end-plate surface by 10° to 30° elevons varied from 85 to 98 percent of the force 
that would be predicted using the 2-D flow analysis of equation (1 5) (refs. 2 
and 32) . 


Heat-Transfer Distributions 

Peak heat-transfer coefficients .- The maximum value of the heat-transf er- 
coefficient ratio, induced by the elevon on the end-plate surface, may be esti- 
mated by using (frcxn ref. 8) 


^peak 

= c (Ml sin 0 - 1) + 0.75 

hu 


where 


1 + 

2 /x\ 
3\6/ 


for 

x < 56 

(16) 

4.3 

V 

+ O.Oli 

.a) 

for 

X > 56 



These linear expressions for c were obtained by simply fairing lines through 
the data presented by Hayes in reference 8. 

An alternate expression for the value of the maximum heat-transfer- 
coefficient ratio, expounded by Scuderi (ref. 41 ) , is 



X 


0.2 


( 17 ) 


^peak 


/Ppeak\°-®^ 

.2 + 0.2 

\ PI / 


1 + - 
L 


where x is measured downstream from where the eleven generated shock wave 
intersects the wing surface and L is the distance from the plate leading 
edge to the eleven hinge line. 

Finally, another expression for the maximum heating amplification that 
is in agreement with data from many sources was presented fairly recently by 
Neumann and Hayes (ref. 5) 


’'peak /PpeakX 


0.8 


I— J 


( 18 ) 


Sample comparisons .- Peak heat-transfer values were calculated using equa- 
tions (1 6) to (1 8) and compared with experimental results for an example case 
(fig. 14). The example case chosen was an end plate adjacent to a 20° elevon 
(refs. 2 and 45) . The elevon hinge line is located 55 cm downstream of the wing 
leading edge (L = 55 cm). The boundary-layer thickness, measured from profile 
schlieren flow photographs when the end plate was not attached, is 0.86 cm at a 
station 7.6 cm downstream of the elevon hinge line (ref. 45). This measured 
boundary-layer thickness agrees with the value calculated by the Anderson 
and Lewis boundary- layer code, using a measured pressure distribution as input 
(ref. 37) . The local Mach number of the flow over the flat-plate wing, upstream 
of the hinge line, is 5.9 (ref. 45) . The peak pressure ratio, calculated using 
equations (13) and (14), is approximately 9.5 (shown in fig. 13). 

These values, used in equations (16), (17), and (18) to calculate pre- 
dicted peak heating yield the following results: 


Equation no.: (16) (17) (18) 

hpeak/'^u‘ 8.76 8.55 6.06 


The equation proposed by Neumann and Hayes (ref. 5) gives results close to 
the experimental values for the example case shown in figure 14. 

Wing-sweep and tip-fin effect s.- Heat- transfer-coefficient distributions 
on the end plate are similar to pressure distributions on the end plate in that 
they are not affected by wing sweep as long as the tioundary layer is turbulent 
upstream of the elevon (refs. 32 and 45) . For low-aspect-ratio elevons, large 
tip fins may result in a larger region of separated flow ahead of the elevon; 
this would result in a forward shift of the shock and the increased-heating 



region on the end-plate surface. For large-aspect-ratio elevens, or for eleven 
deflection angles less than those required to cause separation, tip fins at the 
outboard edge of the elevens have no effect on the heat- transfer distributions 
on the end-plate surface. 

Cylindrical body .- Neumann and Hayes (ref. 5) present extensive interaction- 
flow data obtained for unswept fins mounted on ogive cylinders. They find that 
the location of peak heating is quite close to the fin on the cylindrical sur- 
face, and is similar to the location of peak heating on a flat-plate-fin configu- 
ration (eq. (10)). The magnitude of the peak heating on a cylindrical surface is 
similar to that for a flat-plate-fin configuration, and may be estimated using 
equation (1 8) (ref. 5) . 

In figures 15 and 16, interaction-flow regions on an end plate and a cylin- 
drical body are compared (ref. 45) . The interaction-flow data are for M-] = 5.9 
and an elevon deflection angle of 20°. The extent of the interaction-flow region 
is larger on the cylindrical body than on the planar end plate, and the demarca- 
tion line on the cylindrical-body surface is curved. In agreement with the 
results presented by Neumann and Hayes (ref. 5) , the locations of maximum heat- 
ing are approximately the same for the cylindrical body and the planar end plate. 
The heating magnification is larger for the cylindrical body as a result of the 
undisturbed heat-transfer coefficient on the cylindrical body being somewhat less 
than that on the end plate (ref. 45) . 


CONCLUDING REMARKS 

Pressure and thermal loads induced by trail ing-edge elevens on high-speed 
aircraft surfaces can differ substantially from those that would be anticipated 
using inviscid-flow analyses. Even without flow separation, the control load 
on the elevon is less than that calculated inviscidly. The extent of the 
interaction-flcw region on an adjacent fuselage surface is much larger than pre- 
dictable using inviscid-flow analyses, but the total side force is less. 

The methods presented herein are adequate for preliminary estimates of the 
forces induced on adjacent planar surfaces by deflected elevens. The extent of 
the interaction-flow region, as well as surface pressure distributions, can be 
predicted for planar wing, elevon, and end-plate surfaces. For nonplanar con- 
figurations, such as a cylindrical fuselage adjacent to the elevon in lieu of 
a planar end plate, there are no analytical methods for predicting the extent 
of the 3-D interaction-flow region. 

Heating magnification factors may be estimated, using the methods pre- 
sented herein, as long as the pressure distributions are known. Regions of 
maximum heating, as well as their magnitudes, are estimable using the methods 
given herein. 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
June 26, 1 981 
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Figure 1 Model of hypersonic research airplane. 



Hypersonic research airplane 



Figure 2.- Outlines of typical hypersonic research airplane and 

wing-elevon model. 




(a) No separation upstream of elevon hinge line. 



(b) Separated flow upstream of elevon hinge line. 


Figure 3.- Three-dimensional inter act ion- flow models. 
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Figure 5.- Streamwise pressure-ratio distributions on upswept wing and elevon 
surfaces along elevon midspan line. M-i = 5.9. 
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(a) Cylindrical body and tip fin. e = 20° 


Separation 


Realtachment 


(b) Cylindrical body and tip fin. e = 30°. 

L-81 -165 

Figure 6.- Planform oil-flow pictures for 70° swept wing with various 
attachments and eleven deflection angles. Mi =5.9. 






Reattachment 


(c) Cylindrical body, g 
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Figure 7.- Separation lengths ahead of trailing-edge flaps. Values of 
^sep/'^O increase with decreasing M-] and R for a specific 
eleven deflection angle. 
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Figure 9.- Region of increased pressure and thermal loads on end-plate surface 

for M] =5.9 and e = 20°. 
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Figure 10.- Sketches of pressure and heat-transfer-coefficient distributions 
in interaction- flow region (from ref. 8). 
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Figure 13.- Sample comparison of pressure distributions calculated using 
analytical method (ref. 8) with experimental data (ref. 2) for M-] = 
and e = 20° at a station 7.6 cm downstream of elevon hinge line. 




^ Calculated 

points 

^ Curves fai] 

red through 

>s 

calculated 

points 

0 Experiment 




-7.6 cm 



Peak location 


Shock location 


Hayes (ref. 8) 


• 0 * 0 . \ 


Sender! (ref. 41) 


Neumann & Hayes (ref. 5) 

Separation 

location 


o o o o o 


Peak Shock 


Sep. Onset 


Figure 14.- Sample comparisons of peak heat-transfer coefficients for M-| 

and e = 20°. 
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Figure 16.- Extent of interaction flow and region of high heating on 

cylindrical- body surface. 
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